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The inclusive production of high-p ^ particles (pions) in the beam fragmentation regions 
of high-energy hadronic collisions is driven by the breakup of valence constituents of the 
beam hadrons into their two— body Fock— state components, and their subsequent fragmen- 
tation. We briefly discuss an approach, that allows the consistent inclusion of intrinsic 
and radiatively generated transverse momenta of initial state partons, and describe an 
extension of our approach to nuclear targets. 
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1. Linear k± factorization 

Inclusive particle production at large transverse momenta p± is conventionally 
treated in the collinear factorization approach, well known from textbooks ™ A 
crucial assumption is that p± is the only large dimensionful scale, and the Bjorken- 
x's of colliding partons are not too small. The differential cross section takes a 
factorized form, schematically 

da(pp — > MX) ~ rii eg) rij ® da(ij — > cd) ® D c -,m , (1) 

and is calculable in terms of collinear parton densities u» j, extracted, e.g. from deep 
inelastic scattering, the fragmentation functions as extracted, e.g. from e + e~ annihi- 
lation, and the parton level cross section calculated in pQCD. Here, at leading order 
the parton level final state is a back-to-back dijet system. At higher energies, when 
perturbatively large transverse momenta may satisfy Aqcd < pi < y/s, however 
the collinear factorization is inadequate, and an explicit inclusion of transverse mo- 
mentum degrees of freedom for initial state partons is called for. The pertinent 
formalism - (linear) k± factorization, many applications and some open problems 
are described in Refs. (0 EJ). In the problem of interest here, production of par- 
ticles at large (pseudo)-rapidities in the beam fragmentation region, the relevant 
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target partons will be dominantly gluons, carrying small x (see e.g. Ref. ©). The 
dominant beam partons will clearly be valence quarks, and indeed the parton level 
production process may be seen as a scattering of valence quarks into large p± at a 
Born level, and as a breakup into their quark-gluon component at the level of the 
hrst radiative correction EE1 The distribution of valence quarks is a steep function 
of their momentum fraction — > 1, and the production of high p± forward 
partons must be sensitive to their energy loss in the q — > qg transition. Evaluation of 
the p± -dependent energy loss requires a calculation of virtual radiative corrections 
to the radiationless Born term within the framework of k± factorization, which has 
been accomplished in Ref. (0). We summarize some of its salient features. First, the 
parton level cross section is a linear probe of the target unintegrated gluon den- 
sity, emission of slow gluons is consistently absorbed into the BFKL evolution of 
the target unintegrated glue, second the produced qg dijet system in the breakup 
is not back-to-back, and its azimuthal decorrelation maps the target unintegrated 
glue, third there is a smooth matching, at the leading-log (pj_) -level to the NLO 
of collinear factorization. In Fig. 1 we compare a calculation using a realistic un- 
integrated glue with recent data from the STAR Collaboration^. Notice that here 
no initial state smearing was included, and improvement can be expected, see e.g. 
Ref. (|3J), although especially at lower p± such smearing will depend on a (model 
dependent) extrapolation of the pQCD cross section. A detailed discussion will be 
found elsewhere Q. 

2. Nonlinear k± factorization 

The crucial assumption behind the linear k± factorization discussed above was 
that we could restrict ourselves to the single ^-channel gluon exchange, and could 
neglect the constraints of s-channel partial wave unitarity. Such assumption is war- 
ranted, in a limited energy range, by the finding of only a small fraction of diffrac- 
tive DIS-events at HERA, in the color dipole-proton scattering at HERA energies 
Cei/ctot <C 0.5. Still, unitarity constraints will be looming at LHC energies, even 
in pp collisions, and the fate of k± -factorization in a regime of strong absorption 
is an important issue. It is then particularly fortunate that heavy nuclei provide 
a possible testing ground for unitarity effects. The latter will be enhanced by the 
size/opacity of a nucleus and are controlled by a new large parameter, the satu- 
ration scale Q\ cx A 1 / 3 . A pertinent quantity in the emerging formalism is the 
unintegrated glue of a nucleus, <f>(b, k) (here b is an impact parameter, and k the 
transverse momentum of the gluon) which is calculable as an expansion over multiple 
convolutions of the free-nucleon glue /(k). Important features of the uni nteg rated 
glue for heavy nuclei constructed from a realistic free-nucleon input are a) a 
saturation scale Q\ ~ 0.8+lGeV, b) a large-K Cronin-type antishadowing enhance- 
ment, c) it furnishes a linear k± factorization of inclusive deep inelastic scattering, 
forward single-jets in DIS, and diffractive dijet amplitudes, d) as shown in [HI the 
first steps of its small- X evolution are governed by a certain nonlinear version of 
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BFKL, the Balitsky-Kovchegov equation. While in DIS on the free nucleon target 
azi muth al decorrelations of the forward dijets map out the proton's unintegrated 
gl ue E31 we found in (9) that linear k± -factorization for the dijet spectrum breaks 
completely in the saturation regime. As a result an extension of k± factorization 
to a strongly absorptive nuclear environment has been worked out, and dubbed 
nonlinear k± factorizatio n^^ Space limitations only allow a rough sketch of 
our most relevant results for the induced radiation process of interest: for rapidities 
y > log(l/xA) , xa ~ I/Raw-m (Ra is the nuclear radius, tun the proton mass), 
the transition q — > qg proceeds coherently over the whole nucleus. Multiple scat- 
terings enhanced by the nuclear size/opacity are evaluated in a Glauber-Gribov 
theory with coupled color channels. Nonlinear k± factorization quadratures for all 
the relevant partonic subprocesses, q — > qg, g — > qq, g — > gg have been obtained in 
Refs. (1 1 1 1) - Here we also mention results obtained in the Color Glass Condensate 
approach, which report a similar breaking of linear /cj_ -factorization (for reviews 
and references, see e.g. Refs.{T2$). Many of the nontrivial features derive from its 
color coupled channel aspect. Indeed for the breakup q — > qg two reaction classes 
are of relevance, the first being excitation from lower to higher color multiplet: 
3(§)8d6 + 15. Here the cross section assumes a form 

da{q -> qg(6 + 15)) cx / $ JSJi9 <g> da* ® $fs/, 9 ® $FSi,q , (2) 
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Fig. 1. Invariant cross section for inclusive n° production at = 200 GeV in pp collisions, at 
mean (pseudo-)rapidity r\ = 3.8 of pions. Data from STAR collaboration 6 . The solid and dashed 
lines are for two different sets of fragmentation functions 14 KKP and Kretzer, respectively. 
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which encodes at the same time probablistic p^-broadening through initial and 
final state scattering (the relevant $>fsi,isi being calculable in terms of the free 
nucleon unintegrated glue), as well as coherent distortions of the q — > qg transition, 
which enter the effective cross section da* . The second reaction universality class is 
a rotation within the initial state multiplet, 3®8d3, here we have: 

da(q -» qg(3)) oc (f>(h, K)\ip dist (z, p - k) - ip(z,p- zk)\ 2 . (3) 

There is a prefactor <j>(b, k), but the dependence on <p is highly nonlinear, be- 
cause of the coherently distorted light-cone wavefunction of the q — > qg transition. 
Furthermore here the collinear pole of the q — > qg splitting is manifest in the second 
undistorted wavefunction, and we may interpret it in terms of a nuclear modifica- 
tion of the quark fragmentation. Although it has not been a part of the original 
presentation, we mention, that equations IL'l.'it can be given a simple Reggeon field 
theory interpretation in the spirit of AGK rules, an d can be extended to more refined 
observables, such as topological cross sections EH. These results will be important 
for the calculation of the quenching of forward jets as observed in dA collisions at 
RHIC. 
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